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The decay of fireballs containing the disoriented chiral condensate (DCC) in A+A collisions has 
been analyzed. We found that the high phase-space density and probably the large final fireball 
volume are the favorable factors to extract a DCC signal from the thermal background. Both of 
these factors are expected to take place at RHIC and LHC energies. A coherent pion component 
then can be observed in pion spectra, in Bose-Einstein correlations and in distribution of the ratio 
of neutral to total pions. 

I. INTRODUCTION 

The main goal of the present and future heavy ion experiments is to discover the new states of the matter associated 
with chiral and deconfinement phase transitions. The former may be accompanied by the presence of the disoriented 
chiral condensate (DCC), i.e. correlated region of space-time wherein the quark condensate is chirally rotated from 
its usual orientation in isotopic space (see, for example, the reviews [|l|). The characteristic momentum spread in the 
DCC is fixed by the inverse size of the region. Recent experimental searches of the DCC signals exploit the overall 
distribution of the ratio of neutral to total pions without momentum cut as an experimental matter it is not a simple 
to select low momenta neutral pions. Up to now a clear evidence of the DCC formation has not found. Perhaps, the 
reason is that the DCC signal is strongly contaminated by the thermal background. Therefore the combine strategy 
using the analysis of the one- particle inclusive spectra and correlations function of charged pions on soft momenta 
region could be effective tool to distinguish DCC 

The ability to identify DCC against the background critically depends on the size and phase-space densities of the 
fireballs. For high enough phase-space densities, which can be reached at RHIC and LHC, the "pion laser" phenomenon 
could take place. The conception of "pion laser" was introduced in one of the first papers [|| devoted to the problem 
of symmetrization of pion emission amplitudes at large phase-space densities in finite source. However, the model 
of independent factorized sources (MIFS) [|[-|5) leads to the completely chaotic radiation without coherence. The 
following arguments can explain such a property of the model. The intensive radiation in the MIFS arising at critical 
value of parameters corresponds to the high-temperature Bose-Einstein (BE) condensation of ideal gas in a finite 
system [^. For infinite homogeneous systems, the BE condensate in ideal gas is associated with the phase transition 
where the coherent condensate wave function can be considered as an order parameter. This phase transition is 
conditioned by spontaneous breaking a gauge symmetry. To describe it, one should input a fictitious source into the 
Hamiltonian, go to the thermodynamic limit and after that switch off the source Q. This non-commutative limit 
procedure leads to the ground state with nonzero expectation value of field, i.e., to the coherent state. However, 
for spatially inhomogeneous (effectively finite) systems with the finite number of particles, there is no mechanism 
which results in spontaneous symmetry breaking down. In fact, for finite systems the Hamiltonian of the "pion laser" 
model has to include an evident symmetry breaking term, for instance, an interaction with a classical source. Then 
the phenomenon of "pion laser" might be regarded as an effect of weak coherent signal strengthening in dense boson 
thermal environment. The decay of the DCC is one of the possible mechanisms for classical source to be introduced 
into the Hamiltonian. 

This paper is aimed at finding some experimentally observed consequences of the hypothesis that the DCC can be 
associated with classical sources interacting with quanta of dense bosonic medium. It is worthy to note that though 
the whole picture developed is a simplified version of more realistic scenario, it is believed to include some essential 
features of soft pion production in future experiments at the next generation of heavy ion colliders: RHIC and LHC. 

In Section II the space-time evolution of the DCC field is described by the equations of motion of the linear sigma 
model. Since pions and sigma particles are in a heat bath, the self-consistent Hartree, or mean-field, approximation 
is appropriate for the model description (see, e.g., j|]). According to Ref. ||], the interaction with the surrounding 
heat bath is supposed to result in the finite lifetime of classical pion field, and the corresponding decay width is 
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introduced. It leads to an additional term, describing the interaction of quasipions with classical source, in the 
quasiparticle Hamiltonian. Actually, the DCC exists in a finite hot inhomogeneous system that is surrounded by the 
normal vacuum and hydrodynamically expands with a finite velocity. Thus, it should be stressed that the problem of 
finding solutions of the self-consistent equations is a very complicated one. Here, we use only some general properties 
of such solutions. 

In Section III the quasiequilibrium stage of the evolution is described by the thermal density matrix for inhomoge- 
neous, effectively finite fireball. The DCC decay in thermal bath can violate the chemical equilibrium for quasipions, 
and for this reason the corresponding chemical potentials enters the density matrix. 

When the system loses the local thermal equilibrium, the description based on the thermal density matrix is no 
longer adequate. In Section IV we describe the dynamic system evolution at this freeze-out stage by similar form 



of the Hamiltonian as one at the equilibrium stage (the method is close to 10 1). The main parameters of the 
evolution such as the decay width and the effective time dependent pion masses cannot be, in principle, derived from 
the self-consistent equations at this stage, they are phenomenologic parameters and differ from themselves at the 
quasiequilibrium stage. Although we do not consider here other higher order interactions in the system which could 
result in the BE condensate superfluidity (see, e.g., flit]), we suppose, however, that the freeze-out at high densities 
can be explained on the basis of the superfluidity of the condensate component of pion liquid. 

Results are presented in Section V where signatures of the DCC decay for pion spectra, for HBT correlations and 
NTTo/^TTtot distribution are analyzed. In Section VI we summarize our results and conclusions. 



II. MODEL OF THE DCC DECAY 



The linear sigma model describes the 0(4) chiral fleld (pa = {(7,lf) by means of the Lagrangian with a simple 
nonlinear self-interaction 

C^^d^ip''d''ip'''^^{ip''ip''~v'^)'^ + h<T. (1) 

For the simplicity, let us suppose that in the Hartree approximation there are no cross correlations between different 
components of the field ipa- Then the equations of motion for pion "i" component of the order parameter (classical 
field) ipci = i^i) and for the field fluctuations (quantum quasiparticles) ipq^t can be obtained from the Hartree equation 
for the pion field (pt = tpg^t + ipc,i (see, e.g., ^) 

d^d^'ip.ix) + mUx)ip^{x) = 2A {ip,{x)f , (2) 
here rnf{x) depends on the field averages. Indeed, Eq. (||) can be split into two equations 

d^^d'^ipc.i + ml i(pc,i = ~fc,i, d^d^^ipq^i + mfcpq^i ^ fc,i, ^ = - 2A(y3^ ^, (3) 

where fc.i{x) = corresponds to the mean- field prescription. Eq. (|^) for quasiparticles can be derived from the 
Hamiltonian Hq^t) ~ J cfirHq^iix), 

'HqA^) = + {V^q4x) f + mUx)ipl,ix)) - fcAx)^qAx)^ (4) 

where 7rq,i{x) =Vq,i {x) is the canonically conjugated momentum. 

We suppose that in the mean-field approximation the classical field (pc.i = ^c,i(r, i) — niipc(j,t) is localized in a 
hot and dense region and is a quasistationary state (its energy LUc{t) changes adiabatically) : 

t 

^c(r,t) = iAc(i-,t) +?A,(r,t), ?Ac(r,t) = ^e(r)exp(-z J dt'w,{t')), (5) 

to 

where ni{9, (j)) is the component of the randomly oriented 3- vector n, — 1, we suppose that n does not depend 
on X. The orientation of the vector n can be chosen arbitrarily due to the symmetry of the sigma model Lagrangian. 

The interaction of the classical field with thermal quasiparticles leads to the finite lifetime of the DCC After 
the decay starts, the quantum level is shifted by Sco and the decay width F appears in the perturbation theory 
approximation. Let us neglect Scu and suppose that F(t) = Fi9(t/ — t) + F26(t — tf), Fi ^ mi{x) , F2 = Fi or 
r2 ^ Fi. As it is discussed later, the last prescription corresponds to a rapid disintegration of the system at < > i/ 
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just below the "thermal freeze-out' temperature Tf = T{tf). Then taking the decay of the classical field into account, 
we have: 

t 

=exp(-i J dt'T{t'))^,^,{v,t). (6) 

to 

By using the first equation in (^) and (j^), one can find the deviation from pure mean field prescription in the case of 
DCC decay: 

t 

fcA^) = n,fc{x) = - ^c,^ (x) + exp(-i / dt'T{t')) {x) ^ 0. (7) 



III. THERMAL QUASIEQUILIBRIUM STAGE 

At the initial stage (t^ < t < tf), the system is in a thermal equilibrium and is described by the quasipion 
Hamiltonians Hq,i{t) and the statistical operator p„, 

PnSppn = e-I^.<-W/^, Hg^M = J d\[H,Ax) - ^i^J.^x) ~ ti'{x)J^4x)]. (8) 

Here /i^ are the chemical potentials which are responsible for violation of the chemical equilibrium and J^q.i{x) — 

v'^q^ii^) St y-'q ii^) where (yC^^"* and (/j^ j-* are the positive and negative field components, respectively. The chemical 
potentials appear due to the decay of classical field which accompanies the creation of new quasiparticles with their 
subsequent fast thermalization in the dense environment. To guarantee a finiteness of the system, we include the 
effective phenomenologic "chemical potential" p'{x) in the statistical operator: 

l^'W^A^) = -l^U^Mr), V{r) = m\tf)u.'r' ^ H^i^Lr^ (9) 

Let us suppose that all /i^ are equal to /i and mi{x) — mi{t) ~ m{t). We can disregard the coordinate dependence of 
masses in Eq. (^) because the "trapping potential" V{r) acts as a strong cut-off factor. Note, however, that V{r) is 
not real potential: it is not a part of the Hamiltonian Hq^i{t). At = the quadratic part of ^ in Eq. (||) can be 
diagonalized aXt = tf in the basis ^„(r) = (2^)^n2 (2/)Cn3 (^)j 

4(t/)e„(r) = (-A + V{r) + m\tf))U^), (10) 

with the creation and annihilation operators a\ i{tf), an,i-{tf), where (see, e.g. fl^ ) 

Hm {x^/ m{tin)uj) exp{~x'^m{tf )u;/2) 

tnAx) = — — r , (11 

' ' ((ni!)2"VV™(i/V)^/^ 

elitf) = 2m{tf)u;{ni + ^2 + + ^) + m^{tf). (12) 

At the same time the quadratic part of the quasiparticle Hamiltonian Hq,i{t) can be diagonalized in the momentum 
representation with the operators ai(k,i/) and a\{'k.,tf). The relation between these representations is given by the 
Bogolyubov transformation with the coefficients A„ (k) and i?„ (k) : 

a\{Ktf) )-^\ ^n(k) A„(k) ) \ aUtf) ) ' ^^'^ 

M^)-l(^^+^^]u^), (14) 




Bn{k) - o 1 - 1 en(-k). (15) 
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Eqs. (13) - ( |15| ) can be simply derived if the canonical field variables (coordinate and momentum) are written 
in terms of the creation and annihilation operators in different representations. The coefficients A„(k) and i?„(k) 
satisfy to the condition ^„(A„(ki)j4„(k2) — i3„(ki)i3„(k2)) = 5(ki — k2). The last relation means that it is a 
canonical transformation. It is easy to see that the Fock spaces for aj^ j(t/) and a|(k, t/) are the unitary nonequivalent 
representations of the canonical relations. Physical vacuum is the ground state of the Hamiltonian Hq^i{t) and is 
defined as ai(k, t/) |0) = 0, whereas unphysical vacuum of the "Hamiltonian" iJ^ ^ is defined by a„^i(i/) |0') = 0. So, 
one should renormalizes the ground state of the statistical operator: 

(...)_ = (...)-(0'|...|0'), (...) ^ Sp{...p). (16) 



Below we omit the index " ren" . At large R ^ m{tf) y mjtj) '^^^ calculate the thermal averages of the operators 

a(k,tf), a^(k, i/) in the approximation: yl„(k) 'fn(k), -Bn(k) — > 0. Then the thermal averages at given n in the 
general case fJ. ^ are 

(aj(p,t/)a,(p, i/))^ = Nthip) + I (a.(p, tf))J , (17) 

where 

Ntk{p)^J2 MP'tf))n = ^^J29nUP), (18) 

„ expi,e„ /ij 1 
and gn = / (i^r^„(r)/c(r, t/)/((£„ - ^)^/2e7i). 

IV. FREEZE-OUT STAGE IN THE HEISENBERG REPRESENTATION 

At the final stage (tf < t < tout < oo), the system loses the (local) thermal equilibrium, and the misaligned 
vacuum goes to the normal vacuum with the lifetime r2 oc {tout ~ tf)~^- Here tout is the "physical asymptotic time" 
when aU interactions are neglected and (pc.i(r,tout) — 0, mi(r,tout) — m^, ai{p,tout) = ai^out{p) — Ott, (p, ^out), 
the system evolution is determined by the free pions Hamiltonian. We suppose that the evolution of the system 
during the continuous freeze-out stage tf < t < tout is governed by the Hamiltonian (Q) with the effective masses 
mf{x), as is close to [Q. It should be noted that the coordinate dependence of the effective mass in Eq. (||) ensures a 
finiteness of squeeze-states contributions to spectra. These contributions appear due to the difference in quasiparticles 
masses aX t = tf and t = tout, they are small either in the adiabatic approximation for the mass evolution or when 
mi{r,tf) ~ miir^tout). For the sudden freeze-out scenario we will use the second assumption: mi{r,tf) ~ mi{v,tout). 
In these cases we can make a further simplification supposing mi{x) — mi{t) ~ m{t). Then using the Hamiltonian 
(m, we get the well-known Heisenberg equation's solution for quantum fields which interact with classical source 



ai,outip) ^ aiiP,tf) + riidpitout), (19) 



where 



dp(i) = ^=== / rfi'exp(z / dt"spit"))foip,t'), Bpit) = vV+m^ (20) 
V2ep(t/) J J 

*/ */ 

and /c(p,i) is the Fourier-transformed current fc{x) in (Q). 

In the case of the adiabatic freeze-out F = Fi = r2 ^ 'Ti(t), we have dp{tout) — for all p ii utc < rn. If lOc > m 
then the main contribution has the Breit-Wigner form and reaches its maximal at some momentum p when Sp{t) ~ 
uic{t) at some "average" time point t =i: 



^ n- \ ^ wc(i)V'c(p,t/) / r 

dp{tout)^ ===— — — ^ ^exp(--(t/-to)), (21) 

^2ep{tf)i{ep{t)-uJc{t))- ^ 2 



The sudden freeze-out case corresponds to F2 — » cxd. Then we have 

dpitout) ^ —^{ekVc{p,tf) + ^ <^c (-P,t/)}, (22) 
and in terms of pion field operators, Eq. ( p^ can be rewritten as aT:-{p,tout) — o.-ni{p,tf). 
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V. SIGNATURES OF THE DCC DECAY 



If we relate two stages through Eqs. ( |17D and (|19|), the thermal averages for the operators Ui^outip) at given n are 
written as: 

(4o"t(P)"«.°"*(P))„ = ^*''(P) + |Ko«t(p))„|^ , (23) 



where 



(a»,o«t(p))n = n,C^gn^n{p) + dp{tout))- (24) 



To evaluate inclusive spectra, we should average over all the orientations n. By taking into account that = 
^^e^*"^, no = cos6', the one- particle inclusive spectra for 7r+,7r^,7r° pions are 

If 1 

A^7r+,^-,^o(p) = ^ / A^+,-,o(p)sin6'd6'd(/) = Nth{p) + -A^coh(p), (25) 



where 



iV+,_,o(p) = Nthip) + \n+^-^Q\'^ Ncohip), Ncohip) 



X!-9»^»(P) + d.p{tout) 



(26) 



In order to calculate two-particle inclusive spectra, one has to apply the thermal Wick theorem for fixed n and then 
to average over n. 

To characterize the relative coherent contribution, we consider the ratio " signal to noise" : 

«(P) = %^. (27) 

As it is already mentioned, the classical field is localized in a hot and dense region. The Gaussian length scale of this 
region for rare gas is i?^ ~ R^, while for the dense BE condensate it is ~ — — which is easily seen 

from ^o(r)- Let us suppose for the simplicity that the classical field is localized within the corresponding scale R with 
a Gaussian spatial distribution: 

Mr.tf) = ^r—^ exp(-rV4i?2). (28) 

If the classical field is large enough at the initial moment to and Ti{tf — to) > 1, the DCC decay at the first thermal 
stage can lead to an overpopulation of quasiparticles and hence to appearance of the chemical potential /i 0, as it 
was discussed before. If /z ~ eo, then the average phase-space density is large enough, Upg = {Ncoh + \Nth)/W ^ 1: 
and is mainly determined by the coherent condensate contribution 

nps,coh oc K^— [-r—^ — T ) exp(-ri(t/ - to))- (29) 
£0 \[eo - fJ.)J 

The ratio "signal to noise" can be obtained from (p^, (|2^), (|2 

yieoip) 



Dip) 



Co(p)/(exp(£o-M)-l) 
oc ^£l^ / r^"'P^''^°^ exp(-^(t/ - to)) (X Ky/n^. (30) 



Tf\ eo 2 



The large phase-space density n^g ^ 1 results in the large value 13 (p) ^ 1 due to the weak signal strengthening in an 
overpopulated dense boson medium. This is a "pion laser" phenomenon. This effect is supposed to be experimentally 
observed for pions in the soft energy region: [p^ — m,r) < In rips. In pion transverse energy spectra, this effect 
reveals itself in an essential decrease of the spectrum slope at small energies, , compared to at high energies. 
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Tf. The same effect observes also in a dense non-coherent condensate g-jj], but the thermal freeze-out of such a 
dense system can be hardly explained if there are no coherent attributes (e.g., superfluidity). 

If the overpopulation at the thermal stage does not arise {Ti{tf — to) < 1) and therefore either ~ or it is far 
from the critical value ^ — Eq, the condensate contribution is negligibly small. In the non-relativistic approximation 
and at sufficiently large volume, R^J m[t f)T ^ 1, the thermal part of the radiation can be roughly described (at 



-^('^o ~ l^)\f^J) > 1 ) as the BE gas 



{2Try J exp[(p^/2m(t/) + r^Tf/2R^ + eo - N/Tf\ - 1 

In this case the leading contribution to the coherent emission component Ncohip) is \d-p{tout)\^ , and d-pUout) reaches 
its maximum in the case of sudden freeze-out - which we use here as an illustration. From Eqs. (|2^), ( psj ) and (|3l|), 
we have 

^ \K{tout)? ^ exp(-ri(t, - to)) exp(-p2(2i?2 - (32) 

Some experimental effects can be observed for pions in the soft energy region: (p° — m,^) < (ArriTrR^)^^ which is 
narrower than for "pion laser". Within this region the inverse of energy spectrum slope, (4m^i?^)^^, is considerably 
less than at high energies, Tf. The intcrferometry radius slightly decreases in this gas regime 

The coherence is most directly connected with the intercept of correlation function (CF). For the CF of 7r+7r+-, 
TT~^TT~- pairs the intercepts are 

For -D(p) — > oo one get C~^^{p, 0) = (p, 0) 1.2, and the effective intcrferometry radius squared, with 

C++(p,qe//) - 1 = (C++(p,0) - l)e-''=//«'^//, qlffRlff = 1, (34) 
drastically shrinks at condensation 



"^'/^(P) " l + 2HN^Ap)/NgMy ^^^^ 

Finally, let us discuss the problem of the DCC observation where the distribution of the ratio of neutral to total pions 
/ ~ N^^g/N^r^^^ is analyzed. It is easy to see from Eq. ( p6| ) that the fixed ratio / = N^^^/N^^^^ corresponds to an 
averaging over at a fixed 0. Then using (p6|), we get 

N^og \ + Dcos^e , ^ 

7V,o,, + iV,+ ,e + 7V,-,, l + D ' ' 

where, for example, N^o^e = / N^o^g{p)d^p, iV^o,e(p) = ^/ iVo(p)rf'/>, and D = i/ Ncoh{p)d^p/ J Nth{p)d^P = 
^Ncoh/Nth- Therefore the /-distribution is 



in the interval 1/3(1 + £>)</< (1/3 + D)/{1 + D), and P{f) ^ beyond the above interval. Note, that neither 
thermal nor coherent pion number fluctuations does not taken into account in (^7|) . For a weak DCC signal, D 1, 
the distribution is much more narrow than P{f) = 1/2^/f at Z? — s- oo. Such a weak signal can be observed if the width 
of the distribution (|37|), D/{1 + D), is larger than the width, ^y2j9N^^, of the pure thermal fluctuation distribution 
peaked around / = 1/3. Unfortunately, for rare gas (when the overpopulation at the thermal stage does not arise) one 

can get: D ~ ~ ^ D(0), and so if D(0) is small then the use of f— distributions as signatures 

of the DCC is impeded. However, at large effective sizes of the system which could be created at RHIC and LHC, 
such an analysis of the /—distribution turns out to be effective for a search for a weak DCC signals. Eqs. (E2h, (pi). 
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( p6| ) and ( p^ ) show that for the rare gas case, the value of K^{tf) = K^exp(— ri(t/ — 1^)) characterizes the number 
of coherent pions ^N^oh which, when increasing the volume, is expected to be increased also. If K?{tf) oc B? ■ const, 
then D is not decreasing for larger volume while is decreasing, for example, oc R^^/^ at constant phase-space 

densities. This fact can be used to search for a DCC at RHIC and LHC energies even if the phase-space density is not 
increasing, i.e., if the tendency for the number of pions in the central rapidity area to be approximately proportional 
to the interferometry volume preserves. One can see also that the increase in helps us to distinguish a DCC 

signal over the thermal pion background from spectra and correlations. 

VI. CONCLUSIONS 

Some possible experimental consequences of the DCC decay have been analyzed. Our consideration is based on the 
most general properties of the DCC formation in high energy nucleus-nucleus collisions such as its localization in a hot 
and dense medium, its dissipation by environment and the finite lifetime, the existences of the freeze-out stage where 
the (local) thermal equilibrium is destroyed and the misaligned DCC relaxes quickly enough to the normal vacuum. 
In the picture developed, roughly speaking, a "coherence conservation" low takes place: the large initial classical field 
(DCC) containing virtually a plenty of coherent pions dissipates partially during the thermal stage, but the resulting 
overpopulated medium amplifies the residual weak DCC signal (" pion laser" effect) and makes the number of coherent 
particles to be large again at the freeze-out stage. Such an effect is connected with multiboson effects in the presence 
of classical source and becomes stronger when the freeze-out phase-space density increases. We have found that the 
DCC signals in spectra and correlations can be observed in a narrow energy region of soft pions, the width of this 
region is maximal for the "pion laser" regime. The signal in the distribution of the ratio of neutral to total pions 
can strongly be contaminated by the thermal background. Nevertheless, even if the tendency of constant phase-space 
densities preserves for higher energies, a possibility to detect the DCC signal over the thermal background at so large 
effective volumes as expected in RHIC and LHC heavy ion experiments has been found. 
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